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Both ex situ and in situ spectroscopic ellipsometry (SE) measurements have been employed to 
investigate the effects of HF cleaning on Si surfaces. The hydrogen-terminated (H- terminated) 

Si surface was modeled as an equivalent dielectric layer, and monitored in real time by SE 
measurements. The SE analyses indicate that after a 20-s 9:1 HF dip without rinse, the Si (100) 
surface was passivated by the hydrogen termination and remained chemically stable. Roughness 
of the HF-etched bare Si ( 100) surface was observed, in an ultrahigh vacuum (UHV) chamber, 
and analyzed by the in situ SE. Evidence for desorption of the H-terminated Si surface-layer, 
after being heated to — 550 "C in the UHV chamber, is presented and discussed. This is the first 
use of an ex situ and in situ real-time, nondestructive technique capable of showing state of 
passivation, the rate of reoxidation, and the surface roughness of the H-terminated Si surfaces. 


In the fabrication of ultralarge-scale integrated cir- 
cuits, preparation of native oxide free Si surfaces, the mon- 
itor and control of Si surface passivation and reoxidation 
are extremely important issues. Aqueous HF etching of Si 
surfaces removes the surface oxide and terminates the Si 
surface with atomic hydrogen. 1,2 The hydrogen termina- 
tion retards the Si surface oxidation, and protects the sur- 
face from chemical attack. 3,4 Therefore, the HF cleaning of 
Si surfaces has received increasing attention. However, the 
properties of hydrogen- terminated (H-terminated) Si sur- 
faces under various conditions, and the degree of surface 
passivation and reoxidation are still under investigation. 4 ^ 7 
In this letter, we report results of ex situ and in situ spec- 
troscopic ellipsometry (SE) studies of HF cleaned Si 
(100) surfaces. 

SE is a surface-sensitive, nondestructive optical tech- 
nique used to characterize surface overlayer thicknesses, 
multilayer structures, optical constants of bulk materials, 
and surface changes. 8,9 Ellipsometry determines the com- 
plex ratio of reflectance R p to R S1 defined as 

p=R/R,=\xa{1>)<P, (1) 

where R p and R s are the reflection coefficient of light po- 
larized parallel to (p ) or perpendicular to (j) the plane of 
incidence, and the values of tan(i^) and A are the ampli- 
tude and phase of the complex ratio. 

The pseudodielectric function <e) can be obtained 
from the ellipsometrically measured values of />, assuming 
a two-phase model (ambient/substrate) 8 regardless of the 
possible presence of surface overlayers. For samples with 
surface overlayers or multilayer structures, SE data must 
be numerically fitted according to an assumed model (e.g., 
a three-phase model: ambient/surface overlayer/ 

substrate). Assuming such a model, values of v f,<r (Av |f 4> ; ) 
and A c (Av f ,<l> ; ), defined as in Eq. (1), are calculated. Here 
hv is the photon energy and $ is the external angle of 
incidence. A regression analysis is established to vary the 
model parameters (e.g., layers thicknesses) until the cal- 
culated and experimental values match as closely as possi- 


ble. This process is done by minimizing the mean square 
error (MSE) function as described in Refs. 10 and 11. In 
our study, ex situ and in situ SE measurements were made 
using a Woollam Co. Variable Angle Spectroscopic Ellip- 
someter (VASE™), which was equipped with a beam- 
chopped, rotating analyzer to increase the stray light re- 
jection and signal to noise ratio. 

Si (100) surfaces covered with native oxide from a 
virgin p-type wafer of 14—22 fl cm resistivity were em- 
ployed to study the effects of HF treatments on Si. A piece 
from the wafer was dipped in 9:1 (volume ratio of de- 
ionized water to 49% HF) HF solution for approximately 
20 s with no rinse. SE measurements were made in air, at 
a 75* angle of incidence, before and after the HF dip. Dur- 
ing the measurement, the automated polarizer azimuth an- 
gle P was set to vary with changes of the measured 
^(Av,,<I>y) to minimize experimental errors in the ellip- 



FIG. 1. V and A values of ex situ SE measurements on a Si ( 100) surface: 
(a) before and (b) after the HF treatment The solid line represents the 
experimental data, and the dashed line is the best fit of the SE analysis. 
Assumed models for the SE analysis, in each case, are sketcheJ with the 

plots 
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FIG. 2. Changes in effective thickness of the H-terminated Si surface (H 
surface) as a function of time, after the 20-s 9: 1 HF dip, monitored by the 
ex situ SE measurements. 


someter system. 12 It was set such that ?=$ with P min = 10". 
The spectra] scan, ranging from 3000 to 8000 A, with an 
increment of 100 A, was completed within 2-3 min. The 
SE analysis indicated that an approximately 25 A thick 
native oxide covered the virgin Si surface as shown in Fig. 
1(a). After the HF cleaning, the Si native oxides were 
removed and the Si surface was terminated with H-silicon 
bonds. 1,2 This H-terminated Si surface was modeled as an 
equivalent dielectric layer described by the optical con- 
stants of Si0 2 13 [for the SE analysis shown in Fig. 1(b)]. 
The optica] constants of Si used for the calculations in the 
SE regression analysis are from Ref. 14. A typical “thick- 
ness” of the H-terminated Si surface (H surface) immedi- 
ately after the HF cleaning, as indicated by the SE regres- 
sion analysis, was in the range 14-17 A. Notice that the 
thickness referred to here as an H surface was not the 
actual thickness of the H-surface layer, but the effective 
thickness of the modeled equivalent dielectric layer of 
Si0 2 , which includes possible Si surface microroughness 
after HF etching (as discussed below). The value of this 
thickness as measured by SE was used to monitor the 
changes in the H surface, 

SE measurements were made on this H surface in air at 
room temperature (RT) over a period of several months 
after the HF cleaning. Changes in thickness of the H sur- 
face were monitored as a Function of time as shown in Fig. 
2. The figure shows that the H-terminated Si surface re- 
mained unchanged for over 2 h, and very little reoxidation 
took place within 3-4 days. After two months the reoxi- 
dized Si surface layer saturates at a thickness of —33 A, 
which is thicker than the native oxide before the HF etch- 
ing. The SE study indicates that remarkable surface passi- 
vation has been achieved by the hydrogen termination of Si 
surface dangling bonds, which contributes to the retarda- 
tion of the Si surface oxidation during air exposure. The 
apparent larger thickness of the reoxidized layer provides a 
clue to the Si surface roughness after the HF etching. 

In situ SE was employed to study changes of the 
H-terminated Si surface at elevated temperatures and the 
bare Si surface conditions after HF etching. During the 
in situ measurements, the ellipsometer was attached to an 
UHV chamber, fitted with a pair of low-strain fused-quartz 


FIG. 3. Real time (* 2 ) data from the HF-trcated Si ( 100) surface, inside 
the UHV chamber, at — 550 *C, near the E 2 critical- point energy. 


windows. 11,15 A Si (100) sample was introduced into the 
UHV chamber, immediately after a 20-s 9: 1 HF dip, and 
clamped on a resistor-heater plate that could be rotated 
and tilted by a rotary drive, inside the chamber. Temper- 
atures were measured and controlled by two JT-type ther- 
mocouples, which were calibrated by an infrared (IR) op- 
tical pyrometer. The typical base pressure of the UHV was 
- 1 X 1CT 9 Torr. 

The HF-cleaned Si (100) surface was heated to 
— 550 "C inside the UHV chamber, and real time ellipso- 
metric measurements of one set of ^ and A data were made 
periodically in time (about once every 5 s) at a photon 
energy of ~4.1 eV, corresponding to the critical point en- 
ergy E 2 (i-e., —4.1 eV for Si at -550 "C), 16 while main- 
taining the sample at — 550 *C. These data were converted 
simultaneously to a pseudodielectric function (e) = (ej> 
In general, it has been established 9,11,15 that the 
peak value of the imaginary part (e 2 ) in semiconductors is 
very sensitive to the presence of surface overlayers (eg., H 
surface) and any possible surface defects such as surface 
microscopic roughness. A reduction of the thickness of 
surface overlayer or surface roughness leads to increasing 
values of (€ 2 ). The highest value (e 2 ) corresponds to the 
cleanest and smoothest bulk surface condition. Therefore, 
changes in the H surface induced at elevated temperatures 
can be monitored in real time by measuring the changes of 
<* 2 >- 

Figure 3 shows changes in (e 2 ) in real time at 



FIG. 4. (a) <e 2 > spectrum of HF-treated Si at RT, measured by in situ 
SE, before (dashed line) and after (solid line) being heated to — 350 ”C in 
the UHV chamber. The dotted line is the fit for the rough Si surface 
model, sketched in (b). (c) A sketch of the H-surface structure, contain- 
ing the surface roughness formed by HF etching. 
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— 550°C. An obvious increase of (6 2 ) within the first 
minute at — 550 "C indicates the desorption of 
H- terminated Si surface at — 550 °C. A flat plateau fol- 
lowed the desorption of the H surface indicates a stabilized 
Si surface. The plateau remained unchanged during further 
extended heating. This also suggests that there is no evi- 
dence of the surface quality deterioration (e.g., increasing 
surface roughness induced by heating) at this elevated 
temperature. 15 

SE data for the H-terminated Si surfaces were taken at 
RT before and after the heating, as shown in Fig. 4(a). A 
comparison of the two (e 2 ) spectra at RT, shows a ~7% 
higher <6 2 > peak value at the E 2 critical point energy 
(—4.3 eV at RT) after the heating. This further confirms 
the desorption of the H surface after heating to — 550 "C. 
The evidence for desorption observed in this experiment is 
consistent with results obtained by thermally stimulated 
desorption measurements as described in Ref. 5. 

The RT (€ 2 ) peak value at the E 2 energy after the 
desorption of the H surface, as shown in Fig. 4(a) (solid 
line), is considerably lower than the known value (—46) 
of Si from literature. 14 This indicates a roughened Si sur- 
face induced by the HF etching. The surface roughness was 
modeled as a top Si layer containing 50% voids, as shown 
in Fig. 4(b). The thickness of this rough Si layer was cal- 
culated by the regression analysis, under the assumption of 
the Bruggeman effective-medium approximation 
(EMA). 17 Good fit was obtained with a thickness of 8.5 A, 
as shown in Fig. 4(a) (dotted line). It indicates an approx- 
imate 1-2 monolayer surface microroughness formed from 
the HF etching. 

This model of rough Si surface was applied to charac- 
terize the H surface measured at RT in vacuum, before the 
heating, as shown in Fig. 4(c). In this case, the voids were 
replaced by Si0 2 to model the H termination of the rough 
Si layer, while the thickness of this rough layer was kept 
the same as in Fig. 4(b). On top of the rough Si surface, a 
pure H-surface layer was modeled and described by the 
optical constants of Si0 2 , as we have discussed previously. 
The thickness of —8 A of this pure H-surface layer was 
calculated through the regression analysis. Notice that by 
modeling a rough interface between the top H surface and 
Si substrate, the top pure H surface has an effective thick- 
ness of - 8 A, which is consistent with an expected mono- 
layer H termination of Si surface. 

Same surface-roughness model [as in Fig. 4(c)] was 
used to reanalyze the ex situ RT SE measurements shown 


in Fig. 1. The results show a reduced effective thickness of 
the initial top H surface, after the HF dip, of —5.4 A, and 
a saturated reoxide layer is —21.7 A. This value is quite 
consistent with the thickness of a fresh native oxide surface 
layer. The small difference between the two nitial thick- 
nesses of H surface (i.e., 8.0 and 5.4 A) are likely due to 
different degrees of surface roughness of the two individual 
HF treated Si samples. 

In summary, HF treated Si ( 100) surfaces were inves- 
tigated by ex situ and in situ SE measurements. The SE 
analysis indicated that the Si (100) surface was well pas- 
sivated via a — 20-s 9:1 HF dip with no rinse. Real-time SE 
data showed evidence for desorption of the H-terminated 
Si surfaces at — 550 *C in the UHV chamber. Si (100) 
surface roughness induced by HF etching was observed in 
a UHV chamber, and analyzed by the in situ SE. It was 
shown that a —1-2 monolayer surface roughness was 
formed after our HF etching. This unique surface charac- 
terization, by in situ SE, provides a useful means of study- 
ing and monitoring various Si surface conditions after HF 
cleaning. 
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HIGH FREQUENCY PERFORMANCE OF fr = ( 2 »- O ' 1 

Si, -.Ga./Si, _,Ge,/Si,_,Ge I H BTs /„ - (4*)' R> C t )~ 112 

D. Roscnfeld and S. A. Alterovitz where is given by 


Indexing terms : Bipolar devices , Transistors, Semiconductor 
devices and materials 

The results of a theoretical study of the performance of high 
speed SiGe HBTs is presented. The study includes a group of 
SiGe HBTs in which the Ge concentration in the base is 20% 
higher than that in the emitter and collector (Le. 
y « x + 0*2). It is shown that the composition dependences 
of f T and the are non-monotonic. As the Ge composi- 
tion in the emitter and collector layers is increased, f T and 
fmtx first decrease, then remain constant and finally increase 
to attain their highest values. 


Introduction : The main interest in SiGe alloys stems from 
their potential use in high speed heterostructure transistors for 
microwave and digital switching applications. In the last two 
decades a large number of HBTs based on ITI-V materials 
have been fabricated and characterised. However, it is only 
during the last few years that material growth and processing 
techniques compatible with silicon technology have been 
developed. The new SiGe growth technologies, as well as the 
recently published promising results of Si/Si, _ x Ge^Si HBTs, 
have generated considerable need among material growers 
and device engineers for theoretical estimation of the potential 
performance of these SiGe devices. The composition depen- 
dence of the SiGe HBT figures of merit, such as the current 
gain, the cutoff frequency (f T ) and the maximum frequency of 
oscillation (J has not bccn experimentally or theoretically 
estimated. 

In this Letter, the results of a theoretical study of the com- 
position dependence of f T and / ^ are presented. In the study, 
the high frequency characteristics of a series of npn SiGe 
HBTs with different emitter and collector compositions were 
calculated. To maintain the advantages of an HBT, the Ge 
concentration of the base in each of the HBTs was set to be 
20% higher than that of the collector and emitter composition 
(Le. y = x + 0*2). The composition dependences of the emitter- 
to-collector transit time (t^) and the base-resistance collector- 
capacitance product (R* C c ) were computed and consequently, 
the dependences of f T and f max on the emitter and base com- 
positions were derived. 

Geometry and structure: We followed the design rules present- 
ed in Reference 1, and obtained device dimensions and dopant 
concentrations similar to those obtained in Reference 2. The 
chosen parameters were not meant to project the ultimate 
performance potential of SiGe HBTs and, therefore, the high 
frequency figures of merit are somewhat inferior to those 
recently published [3]. Because the study was comparative in 
nature all devices examined had the same areas (emitter and 
base 2 x 8/im 2 , collector 10 x 8/zm 2 ), same layer thicknesses 
(300, 500, 5000, 5000 A in the emitter, base, collector and sub- 
collector, respectively), same dopant profiles (10 19 , 5 x 10* 
10 17 and 10 2O cm" 3 ) and same operating conditions (IV 
base-collector bias). The HBTs differ in the emitter and collec- 
tor Ge concentrations and in the base concentration, which 
was always set to be 20% higher than the former. 

Device model and material properties: f T and /„ were 
analysed using widely used formulas : 


+ ** + + T<r ( 3 ) 

In eqn. 3 x, and x e are the emitter and collector charging 
times, x* is the base transit time, is the transit time associ- 
ated with the current induced base at high current levels and 
is the transit time through the base-collector depletion 
region (or collector-subcollector depletion region at high 
currents). 

The above transit times depend strongly on the basic 
properties of the materials such as dielectric constants, mobi- 
lities, saturation velocities, bandgaps and dopant concentra- 
tions. The material properties, in turn, depend on the Ge 
concentration and hence it is expected that f T and will 
show a dependence on composition. 

The SiGe material properties and their composition depen- 
dences were adopted from the literature. The resistances, the 
junction capacitances, and the built-in voltages were calcu- 
lated using the expressions given in Reference 4. For the com- 
position dependence of the bandgaps, the intrinsic carrier 
concentrations and the dielectric functions, we use the formu- 
las published in References 5-7, respectively. The bandgap 
discontinuities and their composition dependences were taken 
from Reference 8. The expressions for the carrier mobilities in 
SiGe are based on the theoretical calculation presented in 
Reference 9. For the saturation velocity of electrons in SiGe 
alloys we employed a linear fit between the saturation velo- 
cities of electrons in Si and Ge. The influence of the high 
current on the width and location of the base-collector deple- 
tion region was calculated according to Reference 10. 

Results and discussion: The calculated f T and /***, plotted 
against the collector current density J c for several SiGe HBTs, 
are shown in Figs. 1 and 2. All transistors were identical in 



Fig. 1 Cutoff frequency f T of Si L _ fiejSi , _ fieJSi , -fie x HBTs 
against collector current density , plotted for several values of emitter 
composition 

X t - X , and X h = X, + 0 2 


© Reprinted, with permission, from Electronics Letters, vol. 29, no. 3, 4th February 1993, pp. 260-261. 
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geometry and dopant profile. However, they differed in the 
compositions of the layers. The emitter and collector composi- 
tions were varied from X f = X c =0toX 9 = X € = 0*8 and the 
base composition varied from 0*2 to 1. 



Fig. 2 Maximum oscillation frequency f of Si l _ M GeJSi 1 -fiej 

S/, - x Ge x HBTs against collector current density ; plotted for values of 
emitter composition 

X t = X, and + 0-2 

Figs. 1 and 2 demonstrate the non-monotonic behaviour of 
both f T and / ^ T . For X t ** 0 (pure silicon in the emitter and 
collector), the highest f T and /„ (29 GHz) are obtained at 
collector currents of 72 k A cm and !62kAcm' 2 , respec- 

tively. Increasing the emitter and collector compositions to 0*4 
(and the base composition to 0*6), the highest f T and 
decrease to 23 GHz, and are obtained for collector currents of 
68kAcm" 2 and ISOkAcm -2 . Further increasing the emitter 
composition to 0*8 (and pure Ge in the base), the peak values 



Fig. 3 f T at J e — 72 \ cm 1 and f^ at J e = 162 A cm~ 2 against 
emitter composition 

At these current densities f T and attain their highest values for 
X € - 0*2 

/ T at 72 k A/cm 2 

L,at 162kA/cm 2 


of f T and f^ increase to 62 and 100 GHz, obtained for J c ■» 
72kAcm“ 2 and 162kAcm“ 2 . Fig. 3 emphasises the non- 
monotonic behaviour of f T and by showing the f T and 
values (calculated for J c = 72kA cm“ 2 and J c « 
172kAcm“ a , respectively), plotted against the emitter com- 
position. * 

Conclusion: We have demonstrated the non-monotonic 
behaviour of f T and f As the Ge composition in the 
emitter and collector layers (as well as the base) is increased, 

fr an< * fmax first decrease, then remain almost constant and 
finally increase to attain their highest values. 
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